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In this paper, we use a known duality between expanding and contracting cosmologies to construct 
a dual of the inflationary flow hierarchy applicable to contracting cosnrologies such as Ekpyrotic and 
Cyclic models. We show that the inflationary flow equations are invariant under the duality and 
therefore apply equally well to inflation or to cyclic cosmology. We construct a self-consistent small- 
parameter approximation dual to the slow-roll approximation in inflation, and calculate the power 
spectrum of perturbations in this limit. We also recover the matter-dominated contracting solution 
of Wands, and the recently proposed Adiabatic Ekpyrosis solution. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 

Inflation [TH3] has been the most successful theory in 
explaining the early stages of the universe. It provides a 
simple explanation for why the universe is homogeneous 
and geometrically flat, solving the well-known horizon 
and flatness problems of the standard Big Bang. In addi- 
tion, it provides natural dynamics to produce the primor- 
dial perturbations which are the the source of large-scale 
structure in the universe. The form of the primordial per- 
turbations predicted by inflation arc in good agreement 
with the observations of Cosmic Microwave Background 
(CMB) anisotropics @]-^. Despite its success, there are 
still unanswered questions about inflation, such as the is- 
sue of initial conditions, which have motivated the search 
for alternatives to inflation. The related Ekpyrotic and 
Cyclic scenarios [7HT7] are two of the most widely dis- 
cussed candidates for alternatives to inflation. Both of 
these models require a phase of slow contraction occur- 
ing before a bounce to an expanding phase. Other alter- 
native proposals include String Gas cosmology ^IH HH] , 
which postulates an early cosmological loitering phase, 
and Tachyacoustic cosmology [20l - [22| . which is a model 
with varying speed of sound. 

In inflationary cosmology, perturbations are generated 
inside the Hubble volume at early times. As the universe 
goes through the phase of accelerated expansion, the co- 
moving horizon shrinks and the perturbations stretch to 
super-horizon scales and "freeze out" to become classical 
perturbations. During the subsequent matter/radiation 
dominated era, as the Hubble Horizon grows, the pertur- 
bations fall back into the horizon. In contracting models 
such as Cyclic and Ekpyrotic cosmologies, perturbations 
are produced during a phase of contraction. The sim- 
ilar prediction of the two dynamically different models 
is based on the similar feature of both models, namely 
a shrinking comoving horizon. In inflation a shrinking 
comoving horizon is generated by a period of acceler- 
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ated expansion whereas in cyclic and ekpyrotic models 
it is generated by a phase of quasi-static contraction. 
In both cases, quantum perturbations that were once 
smaller than the horizon size exit the horizon and be- 
come classical. There is in fact a duality relating the in- 
flationary and contracting cases which leaves the comov- 
ing horizon invariant. Therefore, for every inflationary 
model which predicts a nearly scale-invariant spectrum 
of perturbations, there exists a contracting model with 
the same prediction for perturbations. 



This duality between contracting and expanding sce- 
narios suggests that methods of studying the inflationary 
paradigm such as the inflationary flow formalism [531 HI] 
can be applied to contracting scenarios as well. The in- 
flationary flow hierarchy is a convenient framework for 
surveying the parametric space of inflation in a model- 
independent manner. After reviewing flow methods in 
inflation, we derive a dual to the inflationary flow hi- 
erarchy for contracting scenarios and show that the flow 
equations are invariant under the duality relating the two 
cases, and therefore can be employed to study the pre- 
dictions of contracting scenarios in a model-independent 
way. After studying the background evolution we con- 
tinue the discussion by investigating the cosmological 
perturbations in a contracting phase. 



The outline of the paper is as follows: In Sec. [TT] we 
review the known duality between the contracting and 
expanding cosmologies, in Sec. |III[ after a short review of 
the inflationary flow hierarchy, we derive the dual flow hi- 
erarchy for the contracting models and show that the flow 
equations are invariant under this duality. Section [IV| is 
the study of cosmological perturbations in the case of a 
contracting universe. After deriving the general mode 
equations for perturbations, we solve the scalar and ten- 
sor mode equations for two approximate cases of dual 
analogs to power-law inflation and slow-roll inflation for 
the contracting case. We also show that the recently 
proposed "Adiabatic Ekpyrosis" scenario [T7j is a special 
case of family of solutions to the dual flow equations. We 
conclude in section IVl 



II. DUALITY BETWEEN EXPANDING AND 
CONTRACTING COSMOLOGIES 

A number of dualities have been proposed connecting 
expanding and contracting cosmologies. A duality pro- 
posed by Wands links an expanding inflationary universe 
with equation of state w = — 1 to a contracting matter- 
dominated universe, with equation of state w = [25^ . In 
both cases, a scale-invariant spectrum of curvature per- 
turbations is generated. However, such contracting cos- 
mologies are known to be unstable [26, 27j. Another du- 
ality proposed by Brustein et al. [3H] is based on the effec- 
tive action for cosmological perturbations, and likewise 
connects a stable expanding phase to an unstable con- 
tracting phase. In this paper, we concentrate on a third 
duality between expanding and contracting Friedmann- 
Robertson- Walker (FRW) spaces first proposed by Boyle, 
et al. p9t l30] and extended to non-constant equation of 
state by Lidsey |31) . Under this duality, which exchanges 
the scale factor and Hubble parameter, the inflationary 
flow equations remain invariant, which we show in Sec. 
|III| This property makes this duality especially useful, 
since it allows definition of a self-consistent dual analog 
to the slow roll approximation in inflation, and calcula- 
tion of the corresponding perturbation power spectra. 

Consider first the case of an expanding FRW space, 
with zero curvature and metric 



ds^ = dt^ 



\t)d^^ 



(r) (dr^ - dy? 



(1) 



where t is the coordinate time and the conformal time 
is defined by dr = dt/a. The scale factor is increasing 
with time, so Hubble parameter in the expanding case is 
positive and determined by the Friedmann Equation 



H 



1 



3M2' 



(2) 



where Mp = mp\/\/8TT is the reduced Planck mass, and 
p is the energy density of the cosmological fluid. The 
time dependence of the scale factor is determined by the 
equation of state p = wp of the cosmological fluid via the 
Raychaudhuri Equation 



'6ifi^' + ^")^- 



(3) 



In an expanding FRW space, w > —1/3 corresponds to 
decelerating expansion, a < 0, and w < —1/3 corre- 
sponds to accelerating expansion (aka inflation), a > 0. 
From the Friedmann and Raychaudhuri Equations, we 
have 



a dH 
Hlla 



-2(l + «^)- 



(4) 



The size of the cosmological horizon is set by the Hubble 
length H~^, and in comoving units is given by 



(aH)- 



(5) 



From Eq. m, the comoving horizon evolves as 
d{aH)-^ 1 



da 



2a^H 



(1 + 3w) . 



(6) 



Therefore in the case of decelerating expansion, the co- 
moving horizon grows with time. 



d{aH)- 
da 



>Q, w> -1/3, 



(7) 



and in the case of accelerating expansion, the comoving 
horizon shrinks, 



d{aH)- 
da 



<0, w < -1/3. 



(8) 



For a shrinking comoving horizon, quantum modes with 
constant comoving wavenumber k are stretched to super- 
horizon scales and "freeze" as classical perturbations. As 
the equation of state approaches w = —1, the generated 
density perturbations approach a scale-invariant power 
spectrum. 

We wish to construct a contracting cosmology which 
shares with inflation the property that the comoving hori- 
zon size shrinks with time, resulting in the generation of 
a spectrum of primordial superhorizon perturbations. In 
a contracting universe, the scale factor is decreasing with 
time, and the Hubble parameter is negative. 



H = 



Mp/- 



(9) 



In the contracting case, the condition Q is unchanged, 
which suggests a duality between expanding and con- 
tracting universes which exchanges the scale factor a and 
the expansion rate H 

a(0) = H{c^), 

H{<f>) = a((/.), (10) 

so that the comoving horizon ([5| is invariant, 

Uh\ = [aHy^ . (11) 

Here a tilde denotes a quantity in the contracting cosmol- 
ogy. Notice that the dual scenario is itself an expanding 
phase since H{(t)) = a^cf)) > 0, however time reversal 
results in a contracting phase. The equation of state 
transforms under the duality as 



a dH H da 
H da a dH ' 

5 + 9w 



9{l + w)' 



(12) 



(13) 



The comoving horizon in the contracting universe then 
evolves as 



d{aH)-^ 
da 



>Q, w < -1/3, 



(14) 



and 



djaH)- 
da 



<0, w >-l/3. 



(15) 



The contracting dual of inflation (w < —1/3) is then 
w > —1/3. This is the duality of Boyle, et al., who 
showed that for w — const., the primordial density per- 
turbations produced in the contracting dual universe are 
identical to those produced in the corresponding infla- 
tionary universe. (Note that the contracting dual of the 
de Sitter limit of inflation, with w — —1, is the case of 
divergent pressure, w; — )• cjo.) 

In the next section we apply the above duality to gen- 
eral inflationary spacetimes, and show that there is a 
well-defined dual analog of the slow roll hierarchy for in- 
flation. 



III. A FLOW HIERARCHY FOR 
CONTRACTING FRW SPACES 

In this section we define a dual analog for the inflation- 
ary flow equations [23 IMj , which make it possible to in- 
vestigate cosmological dynamics in a model-independent 
manner. After a short review of the method of inflation- 
ary flow we show that the fiow equations are invariant 
under the duality between the contracting and expand- 
ing scenarios and therefore can be applied equally well to 
the case of cyclic cosmology. 

Consider a flat FRW space with energy density domi- 
nated by a single scalar fleld (j) with Lagrangian 



-C = ^g^-'d.cj^d.cl^ -V {<!>). 



(16) 



The equations of motion for the spacetime are given by 
the Friedmann and Raychaudhuri Equations, which for 
a homogeneous field 4> are: 



H^^ 



1 



3A/|, 



3Af|, 



I'P^ + VW 



The equation of motion for the field (j) is: 
4> + 3H,p + V' (0) = 0. 



(17) 



(18) 



If the field evolution is monotonic in time, we can write 
the scale factor a {(p) and Hubble parameter H {(p) as 
functions of the field (f> rather than time. Equations ( 17) 
and ( |18[ ) can then be re-written exactly in the Hamilton- 
Jacobi form 



V{(t)) = MIj,H^{(j)) 
4) = -2MJ,H' {(j)) 



2Ml 






(19) 



We define an infinite hierarchy of "Hubble slow roll pa- 
rameters" [32l [33; by taking successive derivatives of the 
Hubble parameter H with respect to the field (p: 



2Mt 



H'{4>) 



9 H"{(t)) 



Ud 



^ H(cP) ' 

4 H'icP)H"'{cP) 



^\h = {2M 



HH<P) 



2^fM'i0'-'^ 



H{<PY 



(^+l)g(,/>) 



(20) 



where the prime denotes derivatives with respect to scalar 
field (f). We refer to this hierarchy as flow parameters. 
The first fiow parameter e is related to the equation of 
state by: 



(1 



- w) 



a dH 
H~d^' 



(21) 



where we have used Eq. Q. In an expanding universe, 
the comoving horizon shrinks for w < —1/3, or e < 1, 
which is just the case of inflation. For an inflating uni- 
verse, the scale factor increases quasi-exponentially, with 
the Hubble parameter H ~ const, and 



a (X exp 



Hdt 



— „~N 



(22) 



where N is the number of e- folds. The sign convention for 
N above means that N decreases with increasing time, 
and therefore measures the number of e-folds before the 
end of inflation. The number of e-folds N is related to 
the field by 



dN 



-Hdt 



da{(j)) 

1 



V2Mp y^y 



(23) 



where we choose the sign convention that y^ has the 
same sign as H' {(/>), 



-V2Mi 



H'W 



We can then write the parameter e as: 

1 dH 
^^ HdN' 



(24) 



(25) 



Similarly, taking derivatives of the fiow parameters with 
respect to N we can generate an infinite set of differential 
equations relating the parameters: 



de 

dN 



2e(77-e), 



dN =^ ^^'^' 



dN 



[{£~l)v-£e]'XH + '^'+'^XH, (26) 



Taken to infinite order, this hierarchy of flow equations 
completely specifies the evolution of the spacetime [23] : 
once the parameters e (N) , 77 (N) and so forth are known, 
the cosmological dynamics and scalar field potential are 
fixed via the Hamilton- Jacobi Equations (19). 



We wish to construct a similar hierarchy of flow equa- 
tions applicable to contracting universes, utilizing the du- 
ality discussed in Section |ll] Under the duality a ^ H, 
the parameter e transforms as: 



e = 2Mt 



H'{<t>) 



I = 2Mi 



a'{^) 



H{4>)J "VaW 

From the definition of Hubble parameter we have: 

H{(j)) 2M|.a'((/.)' 



(27) 



(28) 



so that the dual parameter e can be seen to be the inverse 
of the fiow parameter e. 



1 
2MJ, 






(29) 



Under a -^ H ^ the dual of the number of e-folds dN 
—da/a is: 






H'{4>) 
1 



V2Mp y?(0)' 



(30) 



where we choose the sign convention that -v/l has the 
same sign as a' {<p) and therefore the opposite sign as 

H'{cj,)/H{<j,), 



V? = +V2Mi 



«'(0) 



1 iJ(0) 



a(0) 2M^H'{4>y 

The number of e-folds and its dual are related by: 

dN = -idN, 
so that we can write e and e as: 

dN 



dN 

'dN' 



dN 



(31) 



(32) 



(33) 



Defining Af as the logarithm of the comoving horizon size. 



A^: 



- In (ai7). 



(34) 



we see that Af is invariant under the duality a ^ H , and 
is related to N and N by: 



dJ\f ={!-€) dN^{l-e) dN, 



(35) 



In an inflating universe, the comoving horizon shrinks 
{dAf/dN > 0) for e < 1. Likewise, in a contracting uni- 
verse, the comoving horizon shrinks for e < 1, so that e 
plays the same role in contracting universes that e plays 
in expanding universes. We deflne the full hierarchy of 
dual fiow parameters by replacing H (0) in Eqs.(20) with 
a(0). 



i = 2Mt 



V 



aiq)) 
2Ml ^^' 



AM 



a(0)' 
^a'{^y%^ 






(36) 



Comparing Eqs. (23) and (30), we see that the dual dN 



has exactly the same form as dN, so that the dual flow 
equations are identical to their counterparts in inflation: 

de „ - , - >N 

2e iv-e). 



dN 
dfj 

dN 



dN 



I' - ev. 



- [{£-l)fj-ie]'~Xa + ^'+'^Xa. 



(37) 



This duality invariance of the flow equations means that 
for every inflationary potential, there exists an equivalent 
dual cyclic potential for which the comoving horizon has 
the same dynamics 29 . 

We can then deflne a dual analog of the slow-roll limit 
of inflation by taking e ^ 1 and fj ^ 1, so that 



de 
— ^ < e, 

dN 



(38) 



and so forth, with higher order flow parameters small. 
The flow parameters can then be expressed approxi- 
mately in terms of the potential by using the Hamilton- 
Jacobi Equations, 



V (</>) = 3M|,ij2 (0) f 1 

= -2MlH'{(t)). 
For the dual-slow roll limit e ^ 1, 

ij2(^) 



1 



3e(0) 



(39) 



v{4>) 



-M 



r4 \TTi / j.m2 _ 1 J, 2 



= -2M|[iJ'(0)] =-; 



(40) 



In the case of slow roll inflation, e ex [H' /H) ^ 1 means 
that the Hubble parameter H and scale factor a evolve 
as: 

1 



H' 



Mil 



V (0) ~ const. 



„Ht 



(41) 



IV. COSMOLOGICAL PERTURBATIONS 



In the dual limit e ex {a' /a) , we have 



A. Scalar power spectra 



H 



const. 
1 

3Mf 



-02 + V (0) 



0. 



(42) 



The dual limit of slow roll inflation is therefore a quasi- 
static contraction, with the scale factor nearly time- 
invariant. Note that the potential in the cyclic dual to 
inflation must be negative. Differentiating Eq. (40) with 
respect to (f> gives 



The scalar perturbations in a flat FRW cosmology with 
a single scalar field can be expressed in terms of a single 
gauge-invariant parameter. Two commonly used param- 
eters are the Bardeen potential <I>, corresponding to the 
potential in Newtonian gauge, and the curvature pertur- 
bation C, corresponding to the curvature perturbation in 
comoving gauge. Assuming for simplicity p = wp, with 
w = const., the gauge- independent variables ( and $ are 
not independent, but are related by: 



y(0) 



H' e' 

~ 2 

H I 

= 2^{l + f,-~e)., 



where we have used the flow equations (37) and 

d 



d 

d(t> 



1 



H' d 

H dN' 



MpV2edN 
Then, to lowest order in the flow parameters. 



1 ^i?(0)V 2 fV{cj,)^' 



(43) 



(44) 



2M|, 



H'icb) 



Ml 



V'{4>) 



(45) 



Using the first flow equation (37), the second flow pa- 



rameter fj can be written in the dual-slow roll limit as: 



1 dl 
2idN 



1 



1 



Ml, V V 



V 



VV" 

{V'f 



(46) 



This differs slightly from the "fast roll" parameter rj de- 
fined in Ref. [TT], but is physically equivalent. We hesi- 
tate to adopt the terminology "fast roll" because the dual 
limit of slow roll infiation is not the limit oi (f)^ ^ V (0), 
but the highly tuned circumstance that the kinetic and 
potential contributions to the energy density almost ex- 
actly cancel. We will use the term dual-slow roll to de- 
scribe this situation. 

In this section, we have defined a dual flow hierarchy 
which is self-consistent to arbitrary order and which ad- 
mits a dual analog of the slow roll approximation, ap- 
plicable to cyclic cosmology. The dual version of the 
quasi-exponential expansion of inflation is a quasi-static 
contraction, with nearly constant scale factor. In both 
cases, the comoving horizon size is shrinking with time. 
The flow equations governing the evolution of the space- 
time are invariant under the duality transformation. In 
the next section, we derive equations for cosmological 
perturbations in the case of a contracting universe. 



c = 



'Sa^il + w) \a'/a^ 



$ 



(47) 



where a prime denotes a derivative with respect to con- 
formal time. Ultimately, the physical variable which is 
most directly related to CMB perturbations is the New- 
tonian potential, i.e. the Bardeen potential $. In the 
case of inflation, it is conventional to express the pertur- 
bations in terms of ( since it is conserved on superhorizon 
scales even when the equation of state is varying, which 
is not true for $. 

Let us consider the relation between the two gauge- 
invariant potentials in more detail. For simplicity we 
again take the case of constant equation of state, w — 
const. For a flat FRW cosmology, the Friedmann and 
Raychaudhuri equations are written as: 



1 



3M|,' 



a 
a 



QMl 



(1 + ?,w)p, 



(48) 



where prime again denotes derivatives with respect to 
conformal time r. Therefore, the scale factor a evolves 



as 



2/3(l+t«) 



(49) 



There are two sets of expanding and contracting solu- 
tions: 



• Expanding: 



> -1/3, 
< -1/3, 



re [0,c5o], 

T e [-cx),0]. 



• Contracting: 



w > -1/3, r e [-00,0], 
w < -1/3, T e [0,CX)]. 



(50) 



(51) 



The equation of motion for the Bardeen potential $ is: 
<^" + 3{l + w){aH)<i>' + wk^^ ^0. (52) 



Using 



we have 



aH 



l + 3u; 



(53) 



$"+^'^^i^r-i$'+7«fc2$ = 0, (54) 

1 + 3w 



which in the long wavelength limit, fc — > 0, reduces to 

(55) 



$"+^ii±^r-i$' = 0. 



l + Sw 
This equation has two solutions 



$0 — const. 



$1 



-(5+3m))/(1+3«)) 



(56) 



Note that in an expanding universe, $1 is always a de- 
caying mode, regardless of w, since 



> -1/3, tG[0,( 



< -1/3, re [-00,0] 



5 + 3w 
l + 3u; 



5 + 3w 
l + 3u; 



<0, 



> 0. (57) 



In the case of inflation, w < —1/3 and r is negative and 
tends towards zero. Therefore <I>i is again a decaying 
mode. Therefore, in an expanding cosmology, the physi- 
cally relevant mode is the constant mode $o- Using Eq. 
(49), the transformation p7| becomes. 



2 [$' + jaH) $] 



(58) 



so that for $ = $0 = const., ( and $ have the simple 
relationship 



5 + 3w ^ 

^ 3(l + u;) 



(59) 



Since ^ is conserved even when w is changing, we can cal- 
culate ^ during inflation, and then calculate the Bardeen 
potential $ j in any later epoch with equation of state Wf 
as El: 



<i>,^ 3(l + ^,) 



(60) 



This is the standard lore for perturbations in inflation. 

The solution $1 behaves quite differently in the case 
of a contracting universe; 



w > -1/3, T e [-00,0] 



w < -1/3, T G [0,CX)] 



5 + 3w 
'l + 3w 



<0, 



5 + 3w 
l + 3w 



> 0. (61) 



Therefore from Eq. (561 $1 is always a growing mode 



irrespective of w and is the mode of physical interest. 



Using Eq. ( 49 1 , we have 



$1 ex 



,3- 



(62) 



From Eq. (47) it is then apparent that $1 does not 
contribute to the curvature perturbation C,- Therefore, 
in a contracting universe, the gauge-invariant variable C 
tracks only the subdominant mode of the Newtonian po- 
tential (^niHnj- In a contracting universe, we must also 
follow the Bardeen potential <i>. It is more convenient to 
work with two new parameters u and v which are related 
to $ and C by [3F|: 



^C 



V = zC, 



where z is given by: 



V2Mpa 



(63) 



(64) 



The variables u and v are not independent, but are re- 
lated by: 



1 /i; 



" = "fc^Sz 



(65) 



(66) 



These transformation equations are equivalent to inde- 
pendent mode equations for u and v. 



Vk = 0, 



k^~ — ]uk^O, 



(67) 



where 9 = 1/z, k is the conformal wave number, and a 
prime denotes derivative with respect to conformal time 
dr = dt/a. We choose to work with a dimensionless 
parameter y instead of r, defined as: 



y 



k 



Writing 



d 



iH 



d 



d7 "dt dN' 

we obtain a differential relation between y and r, 

'1-e 



dy 



kdr. 



so for the second derivative we have 



d^ 


y' <P 


2y{fi~l) d 


dT^ 


T^dy^ 


(1 - if r2 dy 



(68) 



(69) 



(70) 



(71) 



From Eq. (64) and using contracting flow equations 



( 37 1 we obtain^^'^/6' and z" jz in terms of fLow parameters 

2 



z 
z 



wliere 



f ) n-'«' 



aH 



G{ef,,e 



C + W + 6e - 6?77 + i~] - 21 



^ -f)^ ^ 21^ -2fje-fi + 2£ 



(73) 



Substituting (72) and (71 1 into (67) we find 



;a2,.2 



(fvk 



{l-l)y'^-2y{f,-~e) 



dy 



y'i'-F{i,fj,^') 



Vk = 0, (74) 



il^efy' 



dy"^ 



2y (77 - i) 



duk 
dy 



yH^-G 



(2,^,e^) 



uk = 0, (75) 



which are exact mode equations for u and v. In the 
following sections we solve the mode equations in three 
approximate limits: power law contraction, the dual-slow 
roll limit, and Adiabatic Ekpyrosis. In the power-law 
case, we also recover the matter-dominated contracting 
solution of Wands. 



B. Power Law^ Case 

In this section, we consider the dual analog to power- 
law inflation, which has the useful property that it is pos- 
sible to solve exactly for both the background evolution 
and the perturbation. Power-law inflation corresponds to 
taking the first flow parameter e to be exactly constant. 



1 dH 



(76) 



so that the Hubble parameter and scale factor evolve as: 



Hoze'^, aoce"^. 



The comoving horizon then evolves as: 



(77) 



(78) 



which shrinks for e < 1. It is straightforward to show 
that the scale factor evolves as a power law in time. 



a cut 



i/-^ 



(79) 



The dual of power-law inflation for the contracting case 
corresponds to taking the flow parameter e to be exactly 
constant. 



1 da 

e — ^ — const. 

adN 



(80) 



\'^l so that the Hubble parameter and scale factor evolve as: 



cTV 



-N 



a (X e"" , H (X e 

The comoving horizon then evolves as: 



(aH) oc e 



{l-e)N 



(81) 



(82) 



which shrinks for e < 1. Eq. (80) can be simply solved 



for the scale factor as a function of the fleld (p, 

«((/.) (X exp ( J ^-^). (83) 



Likewise, from Eq. (29), we have 



H {cj)) = Ho exp - 



1 



(84) 



/2iMp, 
so that the scalar field potential can be written using the 



Hamilton- Jacobi equation ( 39 ) 



F(0)=(^^)A/2 7j2exp 



?M„ 



(85) 



As in the case of infiation, the scale factor is a power law 



in time, and is the dual of Eq. (79), 

a oc t^. 



(86) 



To express the mode equations ( 74 ) in the power-law 



case, note that for e — const., the fiow equations (37) are 



solved exactly by expressing higher-order flow parameters 
as powers of e. 



e = const., V — ^, C^ = e^i 



^A = e^ (87) 
so the mode equations for v and u are reduced to 

[?V+?(l-2e~)]t-fc = 0, 



{i~~efy' 



d'^Vk 



dy'^ 



{l-ify''^ + [fy'~e]u,^0, 



with solutions 










Vk{y) = Vv 


-<'(A)+«< 


f ey 

vg-i 


uk{y) = Vj/ 


<i.Hi'>(f2L)+<i,HW 


f ey 


where a and /? are given by: 






^^^ 


l-3e 
1-e 


: 





(89) 



1 



(90) 



The integration constants are fixed by specifying appro- 
priate boundary conditions. As in the case of inflation, 
the relevant boundary conditions correspond to canoni- 
cal quantization of the perturbations and the selection of 
a vacuum state, which we take to be the usual Bunch- 
Davies vacuum. For a detailed discussion of the quan- 
tization of the functions v and u, see Ref. |36| : v and 
u are in fact canonical conjugates of one another, and 
the canonically quantized variable is v. Once v has been 
fixed by quantization and vacuum selection, the function 
u (which corresponds to the Bardeen potential $) is de- 
termined by applying the transformation (65). 



The asymptotic form of ( 89 1 for the short wavelength 
limit, ey — > oo, is given by: 



Vkiy) = \/ - i^ie 



— ikr 



C2e' 



(91) 



Selecting the Bunch-Davies vacuum as the vacuum state 
of fluctuations in the ultraviolet limit. 



wj, oc e 



~ikT 



(92) 



requires that C2 = 0. Canonical quantization fixes the 
other constant ci. The canonical commutation relation 
for V is 



[w(x,T),7r(x',T)] =z<53(x-x'). 



(93) 



which corresponds to a Wronskian condition for the mode 

Vk, 

dvt 



*dvk 
The fully normalized solution for v is given by 



My) 




(94) 



(95) 



which agrees with the result of 29 . Having the solution 
for V we use the exact transformation relation between u 



and V (65) to obtain u, 

u{y) = 



1 

fc2 


V — V — 
z 




1 
fc2 


dv dy 
dy dr 


fz'V 

V — 



(96) 



In the power law limit e = ?], so we have 



z 
z 



--C^e-V) 



2/e 
k 



and ( 96 ) reduces to 
Uk{y) = 



1-e 



dv 
dy^ 



(97) 



(98) 



Therefore the fully normalized solution for u is: 



Wfe(j/) = 




(99) 



Note that the argument of the Bessel function is 
iy/ (1 — e) ~ iy. Therefore, mode freezing occurs when 



ey 



Ik 



(100) 



so that in a contracting universe with e <C 1, quantum 
modes freeze out and become classical well inside the 
Hubble length, k ^ (aH). The power spectra for the 
curvature perturbation C and for the Newtonian potential 
$ are obtained by taking the long- wavelength limit. 



Pc(fc) 



2^ 

22/3-4 / p 



M2 



ey->0 

r(3/2); ' 



1-2^ 



,,3-2,3 



(101) 



and 



P<i.(fc) 



2^ 



i2a-2 



4M2 



Hz 

2Mfi 

2 



Uk 



r(«) \ 
r(3/2); 






1-2q 



,l-2a 



(102) 



The spectral index is then calculated by differentiating 
the above equations with respect to k at constant time, 
aH = const. [57] , 



n^- 


- 1 = 


rflnP^ 
dink 


aH—co7ist. 


and 






71$ - 


-1 = 


dlnP,^ 
dink 


aH— const 



== 3 - 2/? = 



l-2a 



-2e 



1 



(103) 



(104) 



In the limit e — > 0, the power spectrum for $ ap- 
proaches scale invariance, while the power-spectrum for 
C approaches the well-known result of a strongly blue 
spectrum n^ = 3. For this case to represent a viable 
cosmological model, we must have mode mixing at the 
bounce, i.e. the growing mode for $ in the contracting 
phase must map to the constant mode in the expanding 
phase. The scale-invariant solution of Wands [5S] is sim- 
ply the limit e = 2/3, which results in n^ = 1, which 
requires no mode mixing at the bounce. 

In the next section, we solve for the perturbations in 
the limit dual to the slow-roll limit of inflation. 



C. Dual-slow roll limit 

The power-law case discussed in the previous section 
has € = fj, and allows exact solution of the background 
and perturbation equations. This is the dual analog of 
the power-law limit of inflation, with e — rj. The slow roll 
limit of inflation is distinct, with e and rj independent, 
but both small. In this section, we consider the dual 
of the slow roll limit, such that the dual parameters are 
independent and small: 



e 7^ 77, e,fj -^1. 



(105) 



Solving the mode equations in the dual-slow roll limit is 
more involved. Unlike the power law case where the term 
proportional to first order derivative of the mode function 
vanishes since e = f^, in the dual-slow roll limit this term 
persists and therefore the Bunch-Davies vacuum is not 
the vacuum state for equation (74 1. For this term to 



vanish, and to obtain a consistent mode equation, we use 
a new parameter x which is related to the dimensionless 
parameter y by: 



dy 



dx 



(106) 



We then get a dual-slow roll mode equation, 

d'^v 
x^{l-2i)-^ + {x^ + 2i-fi)v^0, (112) 



with solution 

Vk{x) > y/x 

where 






\/r^^2i 



C2Hp 



Vl-10g + 4ry 1 
^ 2^r^2? 2 ' 



— X 

(113) 



(114) 



Similar to the power law case, the choice of Bunch-Davies 
vacuum as the vacuum state of fluctuations and canonical 
quantization of v gives the fully normalized solution. 



Vk{y) 




(1) / £2/ 



(115) 



Once again the solution for v can be placed in ( 65 ) to 



find u. The transformation from t; to u in this case is 



■"(2/) 



21 -f] 

? 

ye 



1 — e\ dv 
dy 



(116) 



Therefore there is an identical relation between this new 
dimensionless parameter x and r as that between y and 

T, 



kdr = (1 — e) dy, 
kdr = (1 — e) dx. 

The second derivatives are related by: 

. d^ 



d^ _2l{i)~l) d 
dy^ y (1 — e) dx 



dx^ 



(107) 



(108) 



Hence the mode equation in terms of this new variable 
is: 



(ey)'(l 



e? 



dx^ 



2{iy)[i{i-fj)] 



dv]^ 
dx 



(eyf 



(j,^,!^) 



Vk = 0. 
(109) 



To have a well-formulated Bessel equation we need to 
write the variable ey in the coefficients in terms of x. We 
have a differential relation between the two variables. We 
can expand x as an infinite series in flow parameters by 
integrating by parts. 



ey 



edy = ey 

2~e{l-ff) 
1-e 



dy 

y + 0{e^e^) + 



(110) 



Therefore, to first order in the fiow parameters, 

x~ey. (Ill) 



- = ^ (2e - ,7) 
z ye 



The fully normalized solution for u is then 



Uk{y) = 
where 



fcHl- 



e- 1 



^1 + 6e - Aij 
2^J\-2~e~ 



21- f]. 



(117) 



(118) 



(119) 



The equations for the power spectra for C, and $ are 
the same as in the power law case ( 101|102 ) except that 
the indices a and /3 are different and are given by (114) 
and (119). The corresponding spectral indices are: 



n^ — 1 = 

and 
n$ — 1 = 



dhiR 



dln/c 



= 3-2/3 = 2 + 4g-2?7, (120) 



aH^const. 



dlnP<s> 



dink 



= l-2a^ -4? + 2?7. (121) 

aH^const. 

This can be compared to the result for the case of infla- 
tion. 



n$ — 1 = 7T,(; — 1 = — 4e + 277. 



(122) 



The power spectrum for the Bardeen potential $ is there- 
fore invariant under the duality e — ?> e, 77 — ^ 77 connecting 
inflationary and cyclic cosmology. In the next section, 
we consider the recently proposed "Adiabatic Ekpyrosis" 



scenario. 
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D. Adiabatic Ekpyrosis 

Khoury and Stcinhardt have recently proposed the 
Adiabatic Ekpyrosis scenario |17j . which is a contract- 
ing universe producing a scale-invariant spectrum for the 
curvature perturbation ^. In Adiabatic Ekpyrosis, the 
parameter e evolves with time as: 



_ I dH _ 1 

^^ HdN " t2' 

with a{t) approximately constant, so that 



r 1 1 

z ex ay/e oc - ex — , 

t T 



(123) 



(124) 



which results in a scale-invariant power spectrum for (^. 
In this section we use the dual flow formalism to gener- 
alize this solution to a family of solutions parametrized 
by the second flow parameter ry. 

To construct the required class of solutions, note that 



the flow equations ( 36 ) are solved by the following ansatz: 



i = 'X = 0, 
TJ = const. 



(125) 



This is the dual of the non-slow-roll inflationary solu- 
tions considered in Refs. [571 - HD] . and does not require 
the parameter fj to be small, only constant. This is an 
exact solution for e exactly vanishing, and approximate 
for e ^ 1. We solve for the quantum mode Vk using the 
mode equation (1091, in terms of the variable dx = edy. 
However, for fj = const. ~ C>{1), the expansion (110) 



cannot be truncated after one term as in Eq. ( 111 ). This 



can be seen by writing the first flow equation to lowest 
order in e, 



di 



y 



1 de 2e (f; - e) 



dy l-ldN 



l-l 



2£fj. 



(126) 



The series of partial integrations in Eq. (110 1 is then an 
infinite series in ry. 



(ey)5^(-2,;)" = - 



ey 



2fi 



(127) 



The mode equation ( 109 ) is then 



dx"^ y (1 + 2f]f ^ 
with normalized solution 



Vk =0, 



vu{x)^\MHf\x), 



where 



/3 = 



1 + 4^ 



2(1 + 2fi) 



(128) 



(129) 



(130) 



The transformation (116) in this limit becomes 

'vk {xY 



Uk {x) = - 



V 
1^2r) 



dvk (x) 



dx 



Hi'Hx)., 



where 



1/3-11 



1 



(131) 



(132) 



2(1 + 2^7) 

The spectral indices for C and $ are then, for i) < —1/2, 

2(1 + ^) 
2?7 ' 

and 

2(1 + ^) 



H 



1 = 3-2/3 = 



(133) 



rt$ — 1 = 1 — 2q: = 



(134) 



1-^277 

Therefore for 77 = —1, the power spectra for ( and <& are 
identical and scale-invariant. (This is not in contradic- 
tion to the discussion in Sec. |IV A[ since in Adiabatic 
Ekpyrosis, the equation of state w is changing rapidly.) 
It is straightforward to show that fj = —1 corresponds 
to the case e = 1/e = (l/t^) considered by Khoury and 
Steinhardt, since 



(135) 



H 1 


H ^ 


^ H^ t^ ^^ 


^ l + Ct 


Therefore, for t -> 0, 




d H d 


1 d 


dt I dN 


tdN' 


The first flow equation is then 




de 
dN ^' ' 


-2t^ = -21, 


so that 





f]^l-\^-Y 



(136) 



(137) 



(138) 



and we recover the solution of Khoury and Steinhardt 
[TT] as a special case of the solution (129). For \f]\ ^ 1, 



we recover the dual-slow roll expressions 



n^ 



2(1 + ^) 
\ + 2f] 



2V, 



and 



n$ — 1 



2V 



l + 2fi 



2fi. 



(139) 



(140) 



Adiabatic Ekpyrosis therefore generates a scale-invariant 
power spectrum for cosmological perturbations. How- 
ever, Linde et al. have pointed out that this scenario 
suffers from problems of trans-Planckian curvatures and 
a breakdown of linear perturbation theory [41j . 

In the next section, we consider the generation of ten- 
sor perturbations in contracting universes, and show that 
the tensor spectrum is generically blue, and tensors are 
strongly suppressed relative to the scale-invariant mode 
of the scalar spectrum. 
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E. Tensor perturbations 

In this section, we calculate the spectrum of tensor 
(gravitational wave) perturbations produced in a con- 
tracting universe. To first order, the mode equation for 
tensor perturbations is: 



K + k 



0, 



where 



a 
a 



{any 



21-1 



(141) 



(142) 



Again changing the variable from r to y, the mode equa- 
tion for tensor perturbations is: 



(1-e) 2/ 



dy"^ 



2y (ij - e) 



dy 



[y''e^-e{2~e-l)]hk=0. (143) 



We first consider the case of power law contraction, 
e = i) = const.. In this limit, the mode equation is: 



?\2„,2 



(1-^Ty 



dy'^ 



[y-^l^ - 2i^ -f i] hk = 0, (144) 



and the normalized solution is given by: 

^y \ rr(i) f ^y 



hk{y) 



2\ k \l-l ^ \i-l 



where 



l-3f 



1-e 



(145) 



(146) 



The power spectrum for tensors in the long-wavelength 
limit is then 



^"^Mjli^ 



22t / r(7) 



Ml Vr(3/2) 



hk_ 
a 

e 



1-27 



-\ y3-27. 



The spectral index for tensor perturbations is 



riT 



dlnPr 



din fc 



3-27 



aH— const. 



1 



(147) 



(148) 



The tensor spectrum in a contracting universe is therefore 
strongly blue, with tix ~ 2, and has the same spectral 



index as the curvature perturbation (, ( 103 ), tt-t = n^ — 1. 
This means we can define an analog of the tensor/scalar 
ratio in inflation, r<^, as: 



^C 



Ft _ 16 



(149) 



where we have used Eq. ( 101 1 for the power spectrum of 



the curvature perturbation (^. This is exactly the dual of 
the inflationary expression r = 16e. However, in a col- 
lapsing universe with e ^ 1 , tensors are enhanced relative 
to the curvature perturbation, rather than suppressed as 
in the case of inflation. Nevertheless, both the tensor 
and curvature perturbations are suppressed relative to 
the Bardeen potential $, since 



r$ 



16e 



r(7) 

T{a)J \2{l-\) 



2(a-7) 



y 



2+2Q-27 



Noting that 



4e 



1 



2 (a -7) 

in the nearly scale-invariant limit e ^ 1 , 
r$ ~ IGey^L ,„ — ^ 0. 



y^O 



(150) 



(151) 



(152) 



This makes physical sense, since both Pq and Pt are 
constant modes on superhorizon scales [37], 



^2^3-27 ^ const., 
while P$ is a growing mode, 

i/2yl-2" oc y-2-2-+27. 



(153) 



(154) 



The ratio r^ is therefore constant, but the ratio ?■$ is 
negligible on superhorizon scales, y ^ 1. 

In dual-slow roll limit, the mode equation is given by: 



, d^h 
dx^ 



{l-2~e)x^'^ + {.' + i)h,^0. 



The fully normalized solution is then 



hk{y) 



ey 



1 



£2/ 



e-1 



(155) 



(156) 



which is the same as the solution ( 145 1 in the power law 



case except that the order of Hankcl function 7 is: 



7 



2^jT^Ye 



1 



(157) 



which can be seen as the small-e limit of Eq. ( 146 1. The 
tensor spectral index is then 



•riT 



1 = 



d In Pt 



dlnfc 



= 3 - 27 = 2 + 2g. (158) 



aH^co7ist. 



The tensor/scalar ratio ?■$ is given by (150) with a and 
7 given by (|119) and (157). Having 



2a - 27 = 6? - 2^7, 



(159) 
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in the dual-slow roll limit, e <C 1 and 77 <C 1, the ten- 
sor/scalar ratio reduces to an expression identical to the 
power-law case, 



r$ ~ IGey'^ly^o 



0. 



(160) 



In both the power-law and dual-slow roll cases, the ten- 
sor spectrum has a suppressed amplitude relative to the 
Bardccn potential $, but has constant amplitude relative 
to the curvature perturbation C. The prediction of cyclic 
cosmology for the tensor/scalar ratio is therefore highly 
dependent on the presence or absence of mode mixing in 
the transition from a contracting universe to an expand- 
ing universe. 

Finally, we consider the case of Adiabatic Ekpyro- 
sis. The long-wavelength asymptotic limits for the power 
spectra in (^ and $ are: 



Pc(fc) = 

22/3-4 



M2 



27r2 Z ey^a 

( r(/3) \\ 

\nm)) 



l + 2fi 



1-213 



2ttJ ^ 

(161) 



and 



P$(fc) 



Hz 



2^ 

2^"-^ ( T{a) \ 
AMI 



2M^a 
2 



Uk 



Ti'3/2)J 



iy-^O 

e 
1 + 2r-) 



l-2a 



(162) 



The constants /? and a are given by Eq. ( 130 ) and Eq 



(132 1, respectively. The tensor/scalar ratios for C and $ 
are then 



'^C 



16 
e 



r(7) 
r(/3) 



2 ^~^2(/}-7) 

2 



(1 + 2r)) 



r,\i-2/3 



(1-e) 



~\l-27 



y 



2(/3-7) 



(163) 



and 



r$ 



Pi 



16e 



r(7) 

T{a) 



'^''-''>{l + ^fl)'-'\,-2,+2. 



(1 



-.1-27 



In the scale-invariant limit, e ^ 1 and fj 
sor/scalar ratios are 



■y 

(164) 
-1, the ten- 



re =r$~ 16?y2 1 ^^0^0. (165) 

The tensor/scalar ratio is strongly suppressed for both C 
and $. 



V. CONCLUSIONS 

It has been shown by several authors [25] [29H3T| that 
there exist dualities between the contracting and expand- 
ing cosmologies in their predictions for primordial per- 
turbations. In this paper we have considered the du- 
ality of Boyle et al. and Lidsey [29l |3T], which cor- 
responds to an exchange of the scale factor a{4>) and 
Hubble Parameter H {</>). The comoving Hubble hori- 
zon dn — {a{4>)H[(J3))~'^ is invariant under this duality 
and therefore identical perturbation spectra can be pro- 
duced in inflating and contracting cosmologies. We have 
introduced a self-consistent flow hierarchy for contract- 
ing cosmologies and have shown that the flow equations 
are invariant under this duality. 

After studying the evolution of a contracting back- 
ground in terms of the flow parameters we have in- 
vestigated the evolution of the cosmological perturba- 
tions. The mode equations of perturbations in the gauge- 
invariant Bardeen potential $, corresponding to the New- 
toian potential, are second order differential equations, 
and therefore have two solutions: a constant and grow- 
ing/decaying mode. In the expanding phase the non- 
constant mode is decaying and the physically relevant 
mode is the constant mode. This constant mode of per- 
turbations in the Bardeen potential maps onto the cur- 
vature perturbation C,, so that we can follow C, which is 
conserved on superhorizon scales, to track the perturba- 
tions. However, in a corresponding contracting phase, 
the non-constant mode is a growing mode, which does 
not contribute to perturbations in C,. Therefore, in a con- 
tracting phase we must follow $ as well as C to completely 
describe cosmological perturbations. Having derived the 
general mode equations for perturbations in C, and $, we 
solved the equations in three approximate limits: power- 
law contraction, dual-slow roll, and Adiabatic Ekpyrosis. 
In the power-law limit our result is in agreement with the 
existing results in the literature. In the slow-roll limit, 
our calculation of the spectral indices indicates that the 
spectral index of Bardeen potential is invariant under this 
duality since: 



n$ - 1 = -4e + 277, 



(166) 



which has the same form as the spectral index of the 
Bardeen potential in inflation: 

n$ - 1 = -4e + 2??. (167) 

In contrast, C, has highly blue spectrum, 

nQ-l = 2 + 4i-2f]. (168) 

In the case of Adiabatic Ekpyrosis [T7], both C, and <& 
acquire nearly scale-invariant spectra, 



n^ — 1 = n$ — 1 



2(l + ry) 



l + 2?7 
where the scale-invariant limit is 77 = — 1. 



(169) 
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We also defined the analog of the tensor/scalar ratio r 
in inflation, for contracting models. We have shown that 
the tensor/scalar ratio for the curvature perturbation is 
invariant under the duality, 



''C 



16 

e 



(170) 



However for the Bardeen potential in both the power 
law and dual-slow roll limits, the tensor spectrum is sup- 
pressed relative to Bardeen potential. 



r$ 



where 



ISey^ly^o 



k 
aJI' 



0, 



(171) 



(172) 



In the case of Adiabatic Ekpyrosis, tensors are suppressed 
relative to both ( and $, 



''C = ^* - 16ey^|;^^o 



0. 



(173) 



What we have calulated in this paper are the spectra 
of cosmological perturbations in contracting cosmologies 
before the bounce and subsequent cosmological expan- 
sion. A central issue remains: how do these primordial 
spectra translate into perturbations after the bounce? 
If there is no mode mixing at the bounce, the constant 
mode C in the contracting phase will match to the corre- 
sponding constant mode in the expanding phase, whereas 
the growing mode $ will match to a decaying mode in 



the expanding phase [12H47) . In such a situation, the 
mode of physical interest is (, and the contracting dual 
of slow roll inflation predicts a strongly blue power spec- 
trum, inconsistent with the data. However, the matching 
of modes at the bounce is known to be highly model- 
dependent [48H60] , and the scale-invariant growing mode 
$ can be mapped to the constant mode ^ in the ex- 
panding phase. (It has even been shown that certain 
types of bounce can process the primordial power spec- 
tra in a scale-dependent fashion [6ll-!63l.) Furthermore, 
the predictions of this class of models with respect to the 
primordial tensor/scalar ratio are also dependent on the 
details of mode mixing at the bounce, since the ratio r^ 
is generically strongly suppressed, but the ratio r^ 
can be observably large. Therefore, it is necessary 



171 



170 



to specify in detail the physics of the bounce in order to 
translate the primordial power spectra generated in the 
contracting phase to their counterparts in the expanding 
phase. It remains a tantalizing prospect that bouncing 
cosmologies can provide a fully viable alternative to in- 
flation. 
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